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[i] We have compared six Dst forecast models using 63 intense geomagnetic storms 
( Dst < — 100 nT) that occurred from 1998 to 2006. For comparison, we estimated linear 
correlation coefficients and RMS errors between the observed Dst data and the predicted 
Dst during the geomagnetic storm period as well as the difference of the value of 
minimum Dst (A Dst min ) and the difference in the absolute value of Dst minimum time 
(A t Dst ) between the observed and the predicted. As a result, we found that the model by 
Temerin and Li (2002, 2006) gives the best prediction for all parameters when all 63 events 
are considered. The model gives the average values: the linear correlation coefficient of 
0.94, the RMS error of 14.8 nT, the ADst min of 7.7 nT, and the absolute value of A t Dst 
of 1.5 hour. For further comparison, we classified the storm events into two groups 
according to the magnitude of Dst. We found that the model of Temerin and Lee (2002, 

2006) is better than the other models for the events having — 100 < Dst < —200 nT, and 
three recent models (the model of Wang et al. (2003), the model of Temerin and Li 
(2002, 2006), and the model of Boynton et al. (2011b)) are better than the other three 
models for the events having Dst < —200 nT. 

Citation: Ji, E.-Y., Y.-J. Moon, N. Gopalswamy, and D.-H. Lee (2012), Comparison of Dst forecast models for intense 
geomagnetic storms, J. Geophys. Res., 117, A03209, doi: 10. 1029/201 1JA016872. 


1, Introduction 

[ 2 ] Geomagnetic storms are large disturbances in the 
Earth’s magnetosphere and are produced by enhanced solar 
wind-magnetosphere energy coupling through the magnetic 
reconnection mechanism [Dungey, 1961; Sugiura, 1964; 
Gonzalez et al., 1994]. In general, a geomagnetic storm is 
characterized by a depression in the horizontal (H) compo- 
nent of the geomagnetic field. This depression is mainly 
caused by the ring current encircling the Earth in a westward 
direction and can be monitored by the Dst index [Katnide 
et al., 1998; Daglis et al., 1999; Wang et al., 2003], The Dst 
index is based on measurements from four magnetometers 
near the equator [Sugiura and Kamei, 1991]. It can be cal- 
culated by measuring the horizontal geomagnetic field 
component at four low-latitude ground observatories which 
represent the effects of several current systems in the low- 
latitude geomagnetic field [Akasofu and Chapman, 1972; 
Campbell, 1996; Rostoker et al, 1997; Wang et al., 2003]. 

[ 3 ] Since the Dst is a representative index for indicating a 
geomagnetic storm, there have been several attempts to 
forecast it using real-time solar wind data. Burton et al. 
[1975] found that the Dst can be well modeled using the 
solar wind data obtained from spacecraft as input. Their 
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model predicts the ground-based Dst index solely from the 
velocity and density of the solar wind and the B- component 
of the interplanetary magnetic field (IMF). Following Burton 
et al. [1975], numerous attempts have been made to improve 
the Dst forecast model [Murayama, 1986; Gonzalez et al., 
1994, 2004; Thomsen et al., 1998; Klimas et al., 1998; 
Fenrich and Luhmann, 1998; O'Brien and McPherron, 
2000; Wang et al., 2003; Temerin and Li, 2002, 2006; Xie 
et al., 2008]. Also, the Dst index was predicted by neural 
network methods [Wu and Lundstedt, 1996, 1997 ; Barkhatov 
et al., 2000; Lundstedt et al, 2002; Watanabe et al, 
2003; Pallocchia et al., 2006; Amata et al., 2008] and the 
NARMAX algorithm [Boaghe et al., 2001; Boynton et al, 
2011b]. Some of the above Dst forecast models are con- 
tributing to the prediction of geomagnetic storms via web 
pages as a part of the space weather forecast (e.g., http:// 
sprg.ssl.berkeley.edu/dst_index/ and http://lasp. Colorado. 
edu/space_weather/dsttemerin/dsttemerin.html). 

[ 4 ] The purpose of this study is to compare representative 
Dst forecast models during intense geomagnetic storms 
(Dst < — 100 nT). For this we consider the following six 
models as follows: Burton et al. [1975], Fenrich and 
Luhmann [1998], O’Brien and McPherron [2000], Wang 
et al. [2003], Temerin and Li [2002, 2006], and Boynton 
et al. [2011b]. For the input data for the models, we used 
the solar wind data observed by ACE spacecraft from 1998 
to 2006. For a statistical comparison of the Dst forecast 
models, we estimated linear correlation coefficients and 
RMS (Root Mean Square) errors between the predicted Dst 
data and the observed Dst data during geomagnetic storm 
period. We also obtained the differences between the 
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Days of 2001 


Figure 1. Example of a period of the storm event for the six models. 


observed and predicted values of the following quantities: 
the Dst minimum (A Dst min ) and the time of the Dst mini- 
mum (A t Dst ). Also, we classified the intense geomagnetic 
storm events into four groups depending on the interplane- 
tary structures causing the storms and compared the six Dst 
models for each group. By comparing the statistical models, 
we identify the best Dst forecast model for the prediction of 
intense geomagnetic storms. The paper is organized as fol- 
lows. We briefly explain the data and the six Dst models in 
section 2. We present the results of our study in section 3. A 
brief summary and discussions are given in section 4. 

2. Data and Method 
2.1. Data 

[5] In this study, the observed Dst data were obtained 
from the World Data Center for Geomagnetism, Kyoto 
(http://wdc.kugi.kyoto-u.ac.jp/). The “Kyoto Dst” is calcu- 
lated every hour and is available in real time from this 
website. This real time Dst values are replaced by provi- 
sional values and final values at later dates, which are also 
available at the Web site. Since the time resolution of 
observed Dst data is an hour, we used the ACE magnetic 
field and the solar wind data with the same resolution. The 
ACE spacecraft located in a halo orbit about LI libration 
point ~1.5 million km upstream of the Earth. In this study, 
we considered one hour Dst resolution and did not use 
propagation. The ACE spacecraft data were obtained from 
the Coordinated Data Analysis Web (CD A Web) available at 
NASA’s Goddard Space Flight Center (http://cdaweb.gsfc. 
nasa.gov/cdaweb/). The magnetic field and solar wind 
velocity vectors are given in the Geocentric Solar Magne- 
tospheric (GSM) coordinate system. 


2.2. Dst Forecast Models 

[6] Burton et al. [1975] model (hereinafter referred to as 
the B model) predicts the Dst index based on solar wind data 
obtained from a spacecraft using the following formula: 


where 

Dst* = Dst — b^JPjyn + c. (2) 

[ 7 ] In the above equations Dst* is simply the pressure- 
corrected Dst, from which the contribution of the magneto- 
pause current to Dst has been eliminated by the equation (2). 
When the solar wind dynamic pressure is enhanced, the 
magnetopause moves closer to the Earth, and the currents 
associated with it contaminate Dst. The constant b is a 
measure of the pressure correction and c is a measure of the 
quite-time ring current, magnetopause current, and magne- 
totail current. They determined the constants b and c from 
the linear regression of A Dst and A \fP (b = 16 nT/(nPa) 1/2 
and c = 20 nT). Q(t) is the rate of ring current injection as a 
function of the solar wind electric field, E y = V x B z , where V x 
is an x component of solar wind velocity and fi- is a z 
component of IMF. r is a measure of the ring current decay, 
the value of which corresponds to an e-folding time of 7.7 
hours. Several studies have been performed to improve the 
Dst prediction by either modifying Q or r. Fenrich and 
Luhmann [1998] modified the Burton formula as follows. 
Q(t) was modified to be dependent of the solar wind 
dynamic pressure as well as E v . The value of r corresponds 
to an e-folding time of 3 or 5 hours during the main phase of 
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Table 1 . Start Times, End Times, Minimum Value of Dst, and 
Interplanetary Structure of the 63 Storm Events 


Number 

Start Time (UT) 

End Time (UT) 

Dst m \ n 

Structure 3 

(Type) 

i 

1998/02/17 1200 

1998/02/18 1000 

-100 

nonMC 

2 

1998/03/10 1100 

1998/03/11 1800 

-116 

CIR 

3 

1998/05/04 0200 

1998/05/05 0200 

-205 

SH 

4 

1998/08/26 1000 

1998/08/29 0700 

-155 

nonMC 

5 

1998/09/25 0100 

1998/09/26 0000 

-207 

SH 

6 

1998/10/19 0500 

1998/10/20 0800 

-112 

sMC 

7 

1998/11/09 0300 

1998/11/10 1600 

-142 

sMC 

8 

1998/11/13 0000 

1998/11/15 0400 

-131 

sMC 

9 

1999/01/13 1600 

1999/01/14 2000 

-112 

nonMC 

10 

1999/02/18 0300 

1999/02/19 2100 

-123 

SH 

11 

1999/09/22 2000 

1999/09/23 2300 

-173 

sMC 

12 

1999/10/22 0000 

1999/10/23 1400 

-237 

nonMC 

13 

2000/02/12 0500 

2000/02/13 1500 

-133 

sMC 

14 

2000/04/06 1700 

2000/04/08 0900 

-288 

SH 

15 

2000/05/24 0100 

2000/05/25 2000 

-147 

CIR 

16 

2000/08/10 2000 

2000/08/11 1800 

-106 

nonMC 

17 

2000/08/12 0200 

2000/08/13 1700 

-235 

sMC 

18 

2000/10/13 0200 

2000/10/14 2300 

-107 

sMC 

19 

2000/10/28 2000 

2000/10/29 2000 

-127 

sMC 

20 

2000/11/06 1300 

2000/11/07 1800 

-159 

sMC 

21 

2000/11/28 1800 

2000/1 1/29 2300 

-119 

nonMC 

22 

2001/03/19 1500 

2001/03/21 2300 

-149 

sMC 

23 

2001/03/31 0400 

2001/04/01 2100 

-387 

sMC 

24 

2001/04/11 1600 

2001/04/13 0700 

-271 

sMC 

25 

2001/04/18 0100 

2001/04/18 1300 

-114 

SH 

26 

2001/04/22 0200 

2001/04/23 1500 

-102 

sMC 

27 

2001/08/17 1600 

2001/08/18 1600 

-105 

SH 

28 

2001/09/30 2300 

2001/10/02 0000 

-148 

nonMC 

29 

2001/10/21 1700 

2001/10/24 1100 

-187 

SH 

30 

2001/10/28 0300 

2001/10/29 2200 

-157 

SH 

31 

2002/03/23 1400 

2002/03/25 0500 

-100 

SH 

32 

2002/04/17 1100 

2002/04/19 0200 

-127 

sMC 

33 

2002/04/19 0900 

2002/04/21 0600 

-149 

sMC 

34 

2002/05/11 1000 

2002/05/12 1600 

-110 

SH 

35 

2002/05/23 1200 

2002/05/24 2300 

-109 

SH 

36 

2002/08/01 2300 

2002/08/02 0900 

-102 

SH 

37 

2002/09/04 0100 

2002/09/05 0000 

-109 

CIR 

38 

2002/09/07 1400 

2002/09/08 2000 

-181 

SH 

39 

2002/10/01 0600 

2002/10/03 0800 

-176 

sMC 

40 

2002/10/03 1000 

2002/10/04 1800 

-146 

sMC 

41 

2002/11/20 1600 

2002/1 1/22 0600 

-128 

CIR 

42 

2003/05/29 2000 

2003/05/30 1000 

-144 

SH 

43 

2003/06/17 1900 

2003/06/19 0300 

-141 

SH 

44 

2003/07/11 1500 

2003/07/12 1600 

-105 

CIR 

45 

2003/08/17 1800 

2003/08/19 1100 

-148 

sMC 

46 

2003/11/20 1200 

2003/1 1/22 0000 

-422 

sMC 

47 

2004/01/22 0300 

2004/01/24 0000 

-149 

nonMC 

48 

2004/02/11 1000 

2004/02/12 0000 

-105 

CIR 

49 

2004/04/03 1400 

2004/04/04 0800 

-112 

SH 

50 

2004/07/22 2000 

2004/07/23 2000 

-101 

nonMC 

51 

2004/07/24 2100 

2004/07/26 1700 

-148 

sMC 

52 

2004/07/26 2200 

2004/07/30 0500 

-197 

sMC 

53 

2004/08/30 0500 

2004/08/31 2100 

-126 

sMC 

54 

2004/11/07 2100 

2004/11/08 2100 

-373 

sMC 

55 

2004/11/09 1100 

2004/1 1/1 1 0900 

-289 

sMC 

56 

2004/11/11 2200 

2004/11/13 1300 

-109 

SH 

57 

2005/01/21 1800 

2005/01/23 0500 

-105 

SH 

58 

2005/05/07 2000 

2005/05/09 1000 

-127 

CIR 

59 

2005/05/29 2200 

2005/05/31 0800 

-138 

nonMC 

60 

2005/06/12 1700 

2005/06/13 1900 

-106 

sMC 

61 

2005/08/31 1200 

2005/09/01 1200 

-131 

CIR 

62 

2006/04/13 2000 

2006/04/14 2300 

-111 

sMC 

63 

2006/12/14 2100 

2006/12/16 0300 

-147 

sMC 


Abbreviations are sMC, MC preceded by a fast shock; SH, sheath field; 
CIR, corotating interaction region; nonMC, ICME that does not show the 
signature of a magnetic cloud. 


the storms when E y is greater than 4 mV/m and e-folding 
time 7.7 hours for the remainder of the events. As a result, 
they found that their model (hereinafter referred to as the FL 
model) was much better than the B model. Therefore their 
study implies that during storms the ring current decay is 
indeed bi-exponential and also that the ring current injection 
rate increases during periods of enhanced solar wind 
dynamic pressure. O 'Brien and McPherron [2000] modified 
slightly the Burton formula. The constants b and c were 
determined by a linear fit of the residual phase space offset 
against A \/~P . A linear fit of the residual phase offset after 
the removal of Q to the change of the \/~P yields a pressure 
correction term b of 7.26 nT/(nPa) 1/2 . And they obtained a 
value for c of 1 1 nT. They assumed that the position of ring 
current is controlled by the E y during the ring current 
injection period for southward IMF. Therefore they calcu- 
lated the Q in terms of E y (we refer to this model as the OM 
model). They considered that r also depend on E y . The 
model of Wang et al. [2003] (hereinafter referred to as the W 
model) includes the influences of the solar wind dynamic 
pressure on Q and r. His model finds that the strength of Q is 
proportional to the solar wind dynamic pressure with a 
power of 0.2 during southward IMF and also that r is con- 
trolled by the size of the magnetosphere for northward IMF. 
r decreases as the solar wind dynamic pressure increases, 
interpreted to be the same as the decrease of the magneto- 
spheric size during the periods of northward IMF. Wang 
et al. [2003] compared their model with the FL model, 
OM model, and McPherron and O’Brien [2001]. They 
found that the RMS errors for the predicted negative Dst 
obtained using their model are always smaller than those 
from the other three models. Also, the correlation coeffi- 
cients between the observed and predicted Dst of their model 
are the highest among the four models. Furthermore they 
found that the OM model can be substantially improved, 
especially for intense geomagnetic storms by including the 
influences of the solar wind dynamic pressure. The model of 
Temerin and Li [2002, 2006] (hereinafter referred to as the 
TL model) predicts Dst by defining differently Q and r from 
the other four models. The TL model divides the Q term into 
three parts (dstl , dst2, and dst3), where each Q term is a 
function of the current values of the solar wind magnetic 
field, density, and velocity and of the past values of solar 
wind magnetic field and of past values of dstl and dst2. Q 
terms also have a strong dependence on the angle between 
the solar wind velocity and Earth’s magnetic dipole axis. 
The r terms in the three parts have different time constants 
ranging from 5 days to 1 hour. Temerin and Li compared 
their model with the B model using data for five years 
(1995-1999). As a result, they obtained a RMS error of 
10.6 nT and a prediction efficiency of 67.6%. Using the 
original Burton parameters, they obtained a RMS error of 
1 1.8 nT and a prediction efficiency of 59.7%. The model of 
Boynton et al. [2011b] predicts Dst differently from the other 
models. Their model is based on the NARMAX OLS-ERR 
algorithm that is nonlinear autoregressive moving averages 
model with exogenous inputs (NARMAX) [ Billings et al., 
1989; Boaghe et al., 2001; Balikhin et al., 2001] and an 
orthogonal least squares (OLS) algorithm (hereinafter 
referred to as the NM model). The NARMAX OLS-ERR 
algorithm was employed to identify a single input single 
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Figure 2. CC, RMS, ADst m j n , and I I for the six models. The values in parentheses in Figure 2c are 

the RMS error of A Dst min . 


output model (SISO) for predicting the Dst using a solar 
wind-magnetosphere coupling function deduced by Boynton 
et al. [2011a]. They found the best coupling function for the 
predicting Dst model, C = p 1/2 J' 4 '" B T sm 6 (0/2), where p is 
the dynamic pressure, V is the solar wind velocity, and 9 
is the IMF clock angle ( 9 = tan _1 (5 v /S z )). 

[8] For the input data for the models, we use solar wind 
data from the ACE spacecraft with 1-hour resolution and 
obtain the predicted Dst with the same resolution. For 
comparison, the resolution of the predicted value of the TL 
model was reduced from 5 min to 1 hour by taking the 
average values. We predicted Dst index for 9 years (whole 
period from 1998 to 2006) with first initial Dst value of zero. 
When run the Dst model program for calculating the pre- 
diction value of Dst, we do not use the observed Dst data. 
We use the only previously predicted Dst data. 

2.3. Method of Comparison 

[ 9 ] For comparison we apply the six models to solar wind 
data during 63 intense geomagnetic storms ( Dst< —100 nT) 
from 1998 to 2006 that were identified by Zhang et al. 
[2007] and Echer et al. [2008]. For the quantitative com- 
parisons, we estimated the linear correlation coefficient, the 
RMS error, the A Dst min , and the A t Dst of each event for a 
given model. 

[ 10 ] First, we determined the start times and the end times 
of the 63 storm events using the observed Dst data. Figure 1 
shows the example of period of the storm for the six models. 
The start times are defined as the times when there were 
noticeable decrease of Dst values after the corresponding 
interplanetary shocks (Storm Sudden Commencements, 


SSC) by checking the variation of the solar wind data 
(velocity, density, and magnetic field). The end times are 
defined as the times when the Dst just recovers to above 
—50 nT after the minimum Dst since we are more interested 
in lower values of Dst. The second and the third columns 
of Table 1 show the start times and the end times of the 
63 storm events. 


3. Results 

3.1. Statistical Comparison for All Events 

[ 11 ] We compared statistically the six models by using 63 
intense geomagnetic storm events. Figure 2 shows CC 
(correlation coefficient), RMS, ADst m [ n , and IAfo s ,l aver- 
aged over the 63 storms for the six models. We find that for 
the all parameters, the TL model is better than the other 
models. We can also see that the ADst min values for the W, 
TL, and NM models are much smaller than those of the other 
three models in Figure 2c. Figure 3 shows the event number 
distribution of ADst min for the six models. We see that 
A Dst min values are positive for most of the storms in the 
case of the B and FL model that mostly overestimate the Dst 
minimum values. On the other hand, the event number dis- 
tributions of the W, TL, and NM models are quite symmetric 
and A Dst min most of the events are distributed within 
±50 nT. Therefore, the average values for the W, TL, and 
NM models are small. The RMS errors in A Dst min (shown 
in parentheses in Figure 2c) are also smaller for the W, TL, 
and NM models compared to the other three models. 


4 of 9 


A03209 


JI ET AL.: COMPARISON OF DST FORECAST MODELS 


A03209 



-200 -150 -100 -50 0 50 100 150 200 

A Dst . 



a 

o 

> 

u 

<+H 

o 

u 

0) 

£ 

5 

£ 


-200 -150 -100 -50 0 50 100 150 200 

A Dst . 

min 

40 


30 - 


20 - 



10 - 



-200 -150 -100 -50 0 50 100 150 200 

A Dst . 


40 


30 


20 


10 - 


(d) 


i 1 i 1 — i — 1 — r 


W Model 


n — 1 — r 


j i L 



-200 -150 -100 -50 0 50 100 150 200 

A Dst . 



Figure 3. The event number distribution of A Dst min for the six models: (a) B model, (b) FL model, 
(c) OM model, (d) W model, (e) TL model, and (f) NM model. 


3.2. Dependence on Dst Magnitude 

[ 12 ] For comparison, we classify the storm events into 

following two groups according to the magnitude of Dst: (1) 
Group 1,53 stonn events having — 1 00 < Dst < —200 nT, and 
(2) Group 2, 10 super stonn events having Dst < —200 nT. 
The forth column of the Table 1 shows the Dst minimum 
value of events. 

[ 13 ] Figure 4 shows CC, RMS, ADst mm , and I At Ds , I of the 
two groups for the six models. We see that CC for the TL 
model is the highest among the six models for the group 1 
events, while the CC for the W, TL, and NM models are 
higher compared to the other three models for the group 2 
events (see Figure 4a). RMS values for the TL model are 
smaller than those of the other models for both groups (see 
Figure 4b). The A£>.rf m j n for the TL model is the smallest 


among the six models for the group 1 events, while it is 
the smallest for the W model for the group 2 events (see 
Figure 4c). And, the A Dst m [ n value for the NM model is 
smaller than the other models. Finally, the I At^l value for 
the TL model is the smallest for the group 1 events, while the 
W model has the smallest value for the group 2 events (see 
Figure 4d). From these results, we see that the TL model is 
better than the other models for the stonns with — 100 < Dst 
< —200 nT, while the W, TL, and NM models are better than 
the others for the super storms ( Dst < —200 nT). 

3.3. Dependence on Interplanetary Structures 

[ 14 ] In this section, we classified storm events into four 
groups according to their interplanetary structures causing 
the storms: 27 sMC events (IP shock and MC), 18 SH events 
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Figure 4. CC, RMS, A Dst min , and I Ar /)v ,l of the two groups for the six models. Group 1 includes the 
events having — 100 < Dst < —200 nT and group 2 includes the events having Dst < —200 nT. 


(Sheath Field), 8 C1R events, and 10 nonMC events (non- 
MC type ICME). The interplanetary structures of these 
63 intense geomagnetic storm events were identified by 
Echer et al. [2008]. The fifth column of the Table 1 shows 
the interplanetary structures of events. 

[ 15 ] Table 2 shows the CC for the four groups. We find 
that the CC values for the TL model are the highest for all 
groups. Table 3 shows that the RMS values for the TL model 
are the smallest among the six models for all groups. Table 4 
shows that the ADst m i n values for the TL model are the 
smallest among the six models except for the SH and the 
nonMC group. The ADst m { n for the NM and the OM models 
are smallest for the SH group and the nonMC group, 
respectively. The ADst m] „ for the TL model are all positive 
for the non-MC group events; their average value (about 
16 nT) of the nonMC group events is higher than those (less 
than 10 nT) for the other groups. This fact implies that the 
TL model overestimates Dst for this group. Table 5 shows 
the I Atosi \ of the four groups. We find that I Af/ Jst l for the TL 
model are the smallest among the six models except for the 


Table 2. The CC of the Four Groups for the Six Models 


Model 

sMC 

SH 

CIR 

nonMC 

B model 

0.89 

0.71 

0.77 

0.88 

FL model 

0.83 

0.66 

0.66 

0.81 

OM model 

0.9 

0.8 

0.87 

0.9 

W model 

0.91 

0.81 

0.87 

0.92 

TL model 

0.94 

0.92 

0.95 

0.94 

NM model 

0.88 

0.82 

0.87 

0.89 


nonMC group. The I At£, st \ for the FL model is the smallest 
for the nonMC group. 


4, Summary and Discussions 

[ 16 ] In this paper, we compared six Dst forecast models 
using 63 intense geomagnetic storms ( Dst <—100 nT) that 
occurred from 1998 to 2006. In particular, we examined the 
linear correlation coefficient, the RMS error, A Dst min , and 
A t Ds f The main results from this study can be summarized 
as follows. 

[ 17 ] First, the TL model is better than the other five 
models when all 63 storms are considered for various para- 
meters. Second, the TL model is better than the other models 
for storms with — 100 < Dst < —200 nT, and the three recent 
models (the W, TL, and NM models) are better for the super 
storms {Dst < —200 nT). Third, the TL model is also supe- 
rior to the other models for all types of interplanetary 
structures that cause geomagnetic storms. The detailed 
comparison of these models is summarized in Tables 2-5. 

[is] Comparing statistically the six Dst forecast models, 
we conclude that the TL model is best suited for predicting 


Table 3. The RMS (nT) of the Four Groups for the Six Models 


Model 

sMC 

SH 

CIR 

nonMC 

B model 

51.4 

56.1 

48.2 

39.2 

FL model 

51.7 

50.0 

39.9 

34.2 

OM model 

26.5 

27.6 

19.0 

20.2 

W model 

22.4 

25.4 

16.2 

15.3 

TL model 

18.5 

12.5 

11.8 

11.8 

NM model 

24.1 

26.4 

19.5 

19.6 
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Table 4. 

The ADst min 

(nT) of the Four Groups for the Six Models 

Model 

sMC 

SH 

CIR 

nonMC 

B model 

66.6 

68.0 

69.7 

49.7 

FL model 

63.2 

45.0 

53.7 

52.0 

OM model 

20.3 

19.2 

27.7 

4.6 

W model 

9.6 

9.6 

16.9 

4.8 

TL model 

5.4 

9.0 

2.3 

15.9 

NM model 

14.2 

1.1 

11.7 

5.6 


intense geomagnetic storms. Temerin and Li considered all 
the terms used in the B model and tried to add additional 
solar wind parameters for improvement. The parameters in 
their model are found by minimizing the RMS error between 
the predicted Dst and the observed Dst. In addition, they 
considered several factors different from the other models as 
follows. First, the other models calculate the Dst by using 
the pressure-corrected Dst, the injection term, and the decay 
time. The TL model calculates the Dst by using the six 
terms: dstl, dst2, dst3, pressure term, direct IMF B : term, 
and offset term. The dstl, dst2, and dst3 temis represent 
three current systems within the magnetosphere: the classic 
ring current (trapped ions within the dipole-like field of the 
Earth), the so-called partial ring current, and perhaps the 
magnetotail current, respectively. The physical interpretation 
of dst 3 term is not clear. Second, the four other models use 
only IMF B- components for determining magetospheric 
activity. However, the TL model considers the magnitude of 
the IMF in the B : — B y plane. The IMF B y component also 
makes an effect on the magnetic reconnection and convec- 
tion [Crooker, 1979; Gosling et al., 1985]. Third, while the 
TL model takes into account an annual variation and a 
diurnal variation of Dst, other models do not. The annual 
variation is a function of the angle between the magnetic 
dipole of the Earth and the solar wind velocity and has both 
an annual and a diurnal variation [Cliver et al., 2000; Li 
et al., 2001]. Forth, while the four models use the solar wind 
electric field E y component and/or the solar wind dynamic 
pressure for injection term, the TL model considers the IMF 
clock angle as well as the IMF B- component. The four 
models found the injection functions as g ~ VB, and Q ~ 
VB S P based on the energy-coupling function between the 
solar wind and the magnetosphere of Vasyliunas et al. 
[1982], On the other hand, Perreault and Akasofu [1978] 
and Akasofu [1981] considered that the injection into the 
ring current would occur not only in the southward IMF 
interval but also for a northward IMF. Thus, Q would 
depend on the solar wind-magnetosphere energy coupling 
parameter, e = VB 2 sin 4 (9/2), where B and 9 are the strength 
and the clock angle of the IMF, respectively. However, the 
TL model showed that there is a strong dependence on the 


Table 5. The I AtosA (hour) of the Four Groups for the Six Models 


Model sMC SH CIR nonMC 


B model 

1.8 

2.0 

2.0 

2.7 

FL model 

3.7 

2.0 

5.1 

1.9 

OM model 

1.8 

1.8 

2.2 

3.2 

W model 

1.8 

1.6 

2.3 

2.4 

TL model 

1.5 

1.0 

1.2 

2.3 

NM model 

3.2 

2.1 

2.6 

2.1 


Table 6. CC, RMS, ADst m i n , and Averaged Over the 

139 storms (2003-2010) for the Six Models 


Parameter 

B 

FL 

OM 

W 

TL 

NM 

CC 

0.77 

0.75 

0.80 

0.81 

0.88 

0.77 

RMS 

31.7 

25.0 

16.7 

15.6 

15.0 

16.9 

A Dst min 

41.4 

28.0 

5.8 

5.2 

-3.3 

8.1 

RMS a 

46.0 

35.4 

22.8 

19.7 

17.4 

22.4 

1 ^Dst 1 

2.4 

2.7 

2.5 

2.8 

2.1 

3.2 


a The RMS value of A Dst min . 


direction of IMF such that magnetic activity depends 
approximately on the sixth power of the sine of half the IMF 
clock angle. Balikhin et al. [2010] also showed that the 
coupling functions between solar wind and magnetosphere 
should be proportional to the sixth power of the sine of half 
the IMF clock angle and the NM model predicted Dst by 
using such the coupling function. 

[ 19 ] We have examined the error in percentage to find out 
how well the TL model predicted the Dst minimum value. 
Figure 5 shows the relationship between A.Dst m ,J\Dst min \ 
and IZLsLninl for the 63 storm events for the TL model. We 
find that ^93% (59/63) of the events are located within 
±20%. The RMS error in AZLs7 min /IZAf min l is 0.108. This 
value means that the error in percentage of the predicted the 
Dst minimum value is ~ 1 1 %, which is quite good from the 
viewpoint of space weather forecast. 

[ 20 ] In Figure 5, we examined one noticeable outlier event 
which is the storm that occurred on 20 November, 2003 with 
an error percentage of 30% (printed by arrow in Figure 5). 
This event belongs to the group 2 (see section 3.2). All 
events in the group 2 are located within ±20% except for 
this storm event. As shown in Figure 4, the average value of 
A Dst min for group 2 is 15.3 nT. If we exclude this outlier, 
the average value becomes 3.2 nT. Temerin and Li [2006] 
explained that because they compared the predicted Dst 
(—300 nT) with the provisional version of the “Kyoto Dst” 
index (—472 nT), this event showed a large difference 
between them. Although we compared the predicted Dst 
with the final version of the “Kyoto Dst ” index (—422 nT), 
this event still shows the largest difference among all events. 
On the other hand, the W model very well predicts the Dst 
(—420 nT) of this event; its error is only 2 nT. This event is 
the largest storm of solar cycle 23 and had other peculiarities 
discussed by Gopalswamy et al. [2005]. 

[ 21 ] Temerin and Li found the parameters in the model 
using by minimizing the RMS error between the predicted 
Dst and observed Dst data for eight years (1995-2002). 
They selected optimized parameters for the model by mini- 
mizing the RMS error. Our some events overlap the events 
of Temerin and Li which are used to decide the parameters 
in the model. To avoid some arguments on this overlapping, 
we make another comparison of the six models using 
139 storm events (Dst < —50 nT) from 2003 to 2010. If we 
consider only intensive storms, the number of data is too 
small for statistics. We estimated the linear correlation 
coefficient, the RMS error, A Dst min , and A t Dst . Table 6 
shows CC, RMS, ADst min , and \At Dst \ averaged over the 
all storms for the six models. We can see that for the all 
parameters, the TL model is also better than the other 
models. 
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Figure 5. Relationship between A£>sT min /I.D.sT min l and l/J,v? mm l for the TL model. Here the point symbols 
represent the four interplanetary structures. Super storms ( Dst < —200 nT) are divided by the dotted 
vertical line. 


[ 22 ] Relatively, the W, TL, and NM models well predict 
the super storm events. Among these models, the NM model 
approach possess advantages over that the W and TL 
approaches in that NM provides physically interpretable 
results as demonstrated by Balikhin et al. [2010]. Another 
important point is that the results of the NM model with a 
modified form of the coupling function allow us to revisit 
the analytical derivation of this function and pin-point the 
miscalculation made by Kan and Lee [1979]. 

[ 23 ] Our main purpose is to compare the representative 
Dst forecast models for predicting the intense geomagnetic 
storms and explore suitable models for real-time space 
weather forecast. We conclude that the TL model is the best 
Dst forecast model for most cases, and the W, TL, and NM 
models are the best three models for some super storms 
events. In the near future, we will apply these models to past 
solar wind data in order to check if these models well predict 
the occurrence of geomagnetic storms using contingency 
tables and related statistical parameters [Ji et al., 2010], 
Finally, all these results will be used for real-time Dst 
forecast. 
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